Properties of lightly doped t-J two-leg ladders 
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We have numericaUy investigated the doped t-J ladder us- 
ing exact diagonalization. We have studied both the limit of 
strong inter-chain coupling and isotropic coupling. The lad- 
der scales to the Luther-Emery liquid regime in the strong 
inter-chain coupling limit. In this strong coupling limit there 
is a simple picture of the excitation spectrum that can be 
continued to explain the behavior at isotropic coupling. At 
J = we have indications of a ferromagnetic ground state. 
At a large J/t the ladder is phase separated into holes and 
a Heisenberg ladder. At intermediate coupling the ground 
state shows hole pairing with a modified d-wave symmetry. 
The excitation spectrum separates into a limited number of 
quasiparticles which carry charge -|-|e| and spin i and a triplet 
magnon mode. At half-filling the former vanish but the latter 
evolves continuously into the magnon band of the spin liq- 
uid. At low doping the quasiparticles form a dilute Fermi gas 
with a strong attraction but simultaneously the Fermi wave 
vector, as would be measured in photoemission, is large. The 
dynamical structure factors are calculated and are found to 
be very similar to calculations on 2D clusters. 

PACS numbers: 74.20.Mn, 71.27.4-a, 74.25.Ha 



I. INTRODUCTION 

The properties of strongly correlated electrons confined 
to a ladder (or double chain) and described by t-J or 
Hubbard models haijej-|been the subject of intensive in- 
vestigation rccently.|jl!l The reason lies in the uiHJSiual 
spin liquid nature of the undoped parent system JJw^tll 
Another reason for especial interest is weakly coupled, 
ladders compounds like SrCu203 and (VO)2P207.lla^t£l 
Recent measurements of the magnetic susceptibility and 
the nuclear spin relaxation rate in these materials show 
the existence of a finite spin gap. 

The key question in the current study is the evolu- 
tion of the finite gap in the spin excitation spectrum 
upon doping. The spin gap remains in other spin liq- 
uids systems .Stiid is a sign of strong superconducting 
fiuctuations.EtN 

A recent analysis of the t-J ladder using a mean-field 
theory with Gutzwiller renormalization of the matrix el- 
ements to account for the strong correlations, gave a 
continuous evolution of the spin gap with doping.u The 
short range resonance valence bond (RVB) state evolves 



into a superconductor with modified d-wave symmetry 
within this mean-field approximation. A tendency to- 
wards modified d-wave superconductivity was also found 
in a bosonization approachcl ^and in a recent numerical 
study of the Hubbard ladder .□ 

We have investigated t-J ladders up to a size of 10 x 2 
sites using a Lanczos diagonalization method. First re- 
sults have been published in Ref. ^. Here we report in 
more detail our results for larger lattices including a de- 
tailed investigation of the excitation spectrum, a discus- 
sion of phase separation and the calculation of the super- 
conducting order parameter and of the form factor of the 
Cooper pairs. 

We find clear evidence of hole pairing and a modified 
d-wave RVB state in lightly doped systems in agreement 
with the mean-field theory. An interesting difference 
however is the discontinuous evolution of the excitation 
spectrum upon doping. New "quasiparticle" excitations 
appear carrying both charge and spin. These excitations 
are in addition to a band of magnons which evolve contin- 
uously away from the undoped spin liquid. This separa- 
tion of the excitation spectrum into bound holon-spinon 
quasiparticles and collective magnon excitation contrasts 
with the full spin-charge separation found in a Luttinger 
liquid. 

The t-J ladder Hamiltonian is 

n^-tY^V [c]a,,c,+^.a., + H.c.) V 

-^'E^(4i,-^J-.2,. + H.c.)7' (1) 

j 

where j runs over L rungs, and a (=t, i) and a {— 1, 2) 
are spin and leg indices. The t-J ladder is sketched 
in Fig. |l|. The first two terms are the kinetic energies 
and the J {J') are exchange couplings along the ladder 
(rungs). Unless noted otherwise we set t' — t. The pro- 
jection operator V = Yii a(l ~ '^i.a.T^i.aj) prohibits dou- 
ble occupancy of a site. Periodic or antiperiodic bound- 
ary conditions (PBC, APBC) are used along the ladder. 
The wave vector k = (k^, ky) is consequently well defined, 
kx and ky being the momenta along the ladder and rungs. 
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The transverse momentum ky takes only the values and 
TT, corresponding to bonding and antibonding states. 

At half filling the t-J ladder is equivalent to the 
Heisenberg lajdderp-, which was investigated in earlier 
publicationsJd'BEr EJ The ground state of the Heisenberg 
ladder is,a_slip£.t range RVB state with a spin gap of 
A w J/2iByy at isotropic couphng, J' = J. 

The strong coupling limit J'/ J ^ oo is a good starting 
point to describe the system aa there a simple description 
of the spectrum is available. E£l In that limit, each eigen- 
function of the total system can be written as a direct 
product of one- rung states, which are either spin singlets 
or one of the triplets, and the ground state is that with all 
singlets. The first excited multiplet consists of the states 
with one triplet rung. A small but finite value of J lifts 
the degeneracy of these states. The one-magnon excita- 
tions then form a three-fold spin degenerate band with 
dispersion efe = J' -I- Jcoakx + |(J^/J')(3 — 2cos2fc2;) 
up to second order in J. It has a minimum gap A = 
J' ^ J + \^ at kx — TT.S The momentum perpendic- 
ular to the chains is fcj, = tt. The higher excited states 
form a continuum of excited states and its minimum is at 
k — (0, 0) with energies slightly larger than twice the gap 
2A. With increasing J the collective excitation branch 
crosses into the continuum, but the qualitative descrip- 
tion is still valid. 

In this paper we study the effects of doping holes into 
such a ladder. Although the isotropic case, J'/ J = 1, is of 
most interest, we also study the limit J' 3> J, which can 
be easily understood. In this limit the problem reduces to 
a system of weakly coupled rungs. The properties can be 
continuously followed down to the isotropic point J = J' . 

This paper is organized as follows. In Sec. |l| we briefly 
discuss the occurrence of ferromagnetism in the ladder 
doped with one or two holes at J = and discuss the 
relationship with the occurrence of ferromagnetism in 
two dimensions. Next in Sec. Ill we discuss the pair- 
ing of holes doped into the ladder and the occurrence of 
phase separation. To understand the excitat ion spectra 
we start from the single hole case in Sec. IV and go 
on to the two-hole case in Sec. Section VI discusses 



long range correlations, in particular the interesting ques- 
tion of the symmetry of the pairs and the mapping to a 
Luther-Emery liquid. The single-particle excitations are 
discussed in Sec. |VII| . Over all we find a remarkable sim- 
ilarity between the ladder and 2D clusters. 



II. FERROMAGNETISM FOR J = 



doping in the thermodynamic limit proved to be diffi- 
cult. Actually the ground state of the two-dimensional 
(2D) square-lcittice t-model doped with two holes is not 
ferromagnetic. Lj For finite densities in the thermody- 
namic limit there are contradicting results. Variational 
estimates for the U ^ oo Hubbard model indicate that 
the fully polarized ferroniasnetic state is stable until a 
critical doping Scr- = 0.29.O High temperature series ex- 
pansions by Putikka et al. on the other hand show evi- 
dence that the fully polarized ferromagnetic ground state 
does not survive at any finite doping. Instead they find 
evidence for a partially polarized ferromagnetic state at 
low hole doping. A fully polarized ferromagnetic state at 
finite doping was found only for J < 0.E3 

In this context it is of interest to study the occurrence 
of ferromagnetism in the ladder models. While the proof 
by Nagaokao cannot be applied to the one-dimensional 
chain it is valid for the ladder. The proof relies on the 
existence of closed loops on the lattice. Such loops exist 
in 2D lattices and on ladders, but cannot be formed on 
single chains. The ground state of the ladder doped with 
one hole is thus ferromagnetic. 

We have numerically studied the t-ladder with L = 
2, 3, . . . , 10 rungs, doped with two holes. In Fig. || we 
show the ground state energies of the ladders for both 
PBC's and APBC's. We find that the ground state is 
always ferromagnetic for APBC and an even number of 
rungs and for PBC and an odd number of rungs. For 
the other boundary conditions the ground state is a spin 
singlet. 

An important point is that the ferromagnetic state al- 
ways has the lower energy for a ladder with at least four 
rungs. The singlet state is very close in energy and de- 
serves a more detailed investigation. In Fig. ^ we plot 
the real-space spin correlations (S'^(O)S'^(r)) of the low- 
est singlet state of the L = 10 ladder. These spin cor- 
relations show that the singlet state actually consists of 
two ferromagnetic domains with opposite magnetization. 

The results show clear evidence for a ferromagnetic 
ground state of the i-ladder (L > 4) doped with two 
holes. In the thermodynamic limit however two holes is 
not a finite density. Extrapolations of our small-cluster 
results at finite doping to the thermodynamic limit L — > 
oo are hard to obtain. But one may speculate that the 
existence of a ferromagnetic ground state of the t-ladder 
with 2 holes and L > 4 could indicate a ferromagnetic 
state for dopings 5 < 6^ ~ 0.25. SimilaE-J-esults were 
obtained by Hirsch and Miiller-Hartmann.a 



The t-model {t-J model with J = 0) is equivalent 
to the infinite-// Hubbard model. In single chains the 
ground state of the i-model is degenerate in the spin de- 
grees of freedom. In two dimensions on the other hand 
the ground state of the t-model doped with one hole is 
ferromagnetic. cj This is called the Nagaoka effect. 

The extension of the proof by Nagaoka to finite hole 



III. HOLE PAIRING AND PHASE SEPARATION 

A. Hole Pairing 

In this section we will discuss the pairing of holes doped 
into the t-J ladder and the occurrence of phase separation 
at large values of J/t. We will start from the simple limit 
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J' ^ J,t. In this Hmit the undoped ladder consists of 
weakly coupled rungs, as is sketched in Fig. §(a). 

In this limit two holes doped into the ladder will go 
onto the same rung in order to minimize the number of 
broken singlet bonds. This state is graphically shown in 
Fig. 11(c) . In order to study the occurrence of hole pairing 
at smaller values of J' and down to the isotropic point 
J = J' we calculate the binding energy and the hole- 
hole correlation function. We find that even at isotropic 
coupling the holes still form a bound pair, although the 
pair is more spread out there. 

The binding energy Eb is defined as 



2E, 



G.S. 



(2L - 1) - SG.s.(2i) - Sg.s.(2L - 2), (2) 



where Eq,s.{N) is the ground state energy for N elec- 
trons, the boundary conditions are chosen between PBC 
and APBC to give the lowest energy. 

In the large J' region the binding energy can easily 
be estimated. A single hole doped into a Heisenberg 
ladder breaks one bond with energy loss J', but can 
gain kinetic energy —t along the ladder (see the next 
section for details) and —t' along the rung. It follows 
that Eg,s.{2L - 1) « Eg,s.{2L) + J'-t-t'. Two holes 
on the same rung also break one bond, but the kinetic 
energy of such a bound pair is much smaller, of order 
—At^/J', as will be calculated later. Thus we estimate 
i?G.s.(2i — 2) « i?G.s.(2i) + J', and a binding energy: 



Eb^ J' -2t- 2t' for J' > J, i, t' 



(3) 



Figure ^ shows as a function of J'. It remains 
positive and thus shows binding down to the isotropic 
value, Jjt — J' /t = 0.3. The same holds for a larger 
J/t = 0.5. 

Additional evidence for pairing is provided by the hole- 
hole correlation functions 



(nh(O)nh(r)) = ((1 - n,j_a)(l - rii+r.a')), 



(4) 



measured on the same leg a — a' and on different legs 
a ^ a' in the ground state. They are plotted in Fig. g 
for J/t = 0.3 and J'/J = 1 and 10. For J'/J > 1 the 
two holes are predominantly on the same rung and the 
correlation function shows a clear exponential decay. At 
the isotropic point the pair is more extended. The max- 
imum of the correlation function is now at a distance 1 
along the legs and on different legs, but it again decays at 
large distances. We can calculate the size of the hole pair 
by fitting the inter-chain correlations to an exponential 
form (7ih(0)7T.h(r')) ^ e^''/'^ -t-e"^^^'")/^ for the two largest 
distances, L/2 and L/2 — 1. The inset of Fig. ^ shows 
the size ^ of the hole pair as a function of the inter-chain 
couphng J' /t. The pair is very tightly bound for J' ^ J. 
At the isotropic point the pair is still bound, with a diam- 
eter of about two lattice spacings. Note the oscillation of 
the radius with respect to L. The size seems to converge 
to a value in between the L = 8 and the L = 10 result at 
the isotropic point. 



B. Effective boson model for the large J' limit 

We may say that the system belongs to the Luiher- 
Emery universality class of ID correlated systemsjHj in 
the sense that the spin excitations acquire a finite gap 
while the charge excitations remain gapless. In the limit 
of large J', however, the picture that tightly bound hole 
pairs are moving in a background of singlet rungs is more 
appropriate than weak coupling approaches like g-ology. 
Considering these hole pairs as hard core bosons, we can 
determine the long-range correlations by a mapping to 
an effective boson model. 

The pair hopping matrix element to second order in 
perturbation theory is 



t* = 



2t' 



J' -At'^/J'' 



(5) 



There is a weak attraction V* between two hole pairs on 
neighboring rungs, which again to second order takes the 
form 



V* 



J 
2 



4^2 



3J^ 

8 J' ^ J'-4t'2/J'' 



(6) 



where the first, attractive, term comes from the charge 
part of the J-term in the Hamiltonian. As t* ,V* ^ J' 
we can map the low-energy part of the t-J ladder onto an 
effective hard-core boson model on a chain with nearest 
neighbor interaction: 



H* 



-t 



J2{b}b,+i + H. c.) +V*J2 N^N, 



(7) 



where the hard-core boson creation operator Bj creates 
a hole pair at the rung i and Ni = BjBi is its number 
operator. There is a hard-core repulsion since only one 
hole pair can be created on any given rung. 

Our effective boson model is equivalent to the XXZ- 
model in a magnetic field, which has been solved exactij4 
by a bosonization approach and conformal field theory.E£l 
For V* < —2\t*\ the system is phase separated. This is 



the case for J' > Jpg, where 



PS 



16 r 



2+«U 



(8) 



again to second order perturbation theory. For physically 
reasonable values of J/t phase separation occurs only at 
very large values of J': Jpg/t = 53.2 for J/t = 0.3 and 
Jpg/t = 31.8 for J/t = 0.5. Note that the dominant 
attractive part of the interaction comes from the charge 
part —jJrij^anj+i^a of the J-term. 

Next we will discuss the region where the system is not 
yet phase separated but J' is still large {J,t J' < Jpg). 
There we can determine the dominant correlations from 
the effective boson model. The correlation exponents 
have been calculated indirectly by Bethe ansatz.Ej Both 
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the charge density wave correlations and the supercon- 
ducting correlations show a power- law decay at large dis- 
tances: 



(NrNo)^ const. X r + const, x cos(2fci?r) r 



(9a) 
(9b) 



The superconducting correlations {BIBq) are dominant 
if Kp > 1 . This is the case for most of the phase diagram, 
except for the phase separation regime at V* < —2t*. At 
quarter filling p — 1/2 and for V* > 2t* the system is 
in the Ising-limit and shows a long range charge density 
wave ground state. At fillings close to that line andjJpr 
large V* > 2t* there is a small region where Kp < I.Ej 

In our_effective model we have V* < Q and there 
Kp > 2.l£| We are thus always in the region of domi- 
nant superconducting correlations. Even neglecting the 
attractive charge part of the J-term we are still in the 
superconducting regime where Kp > 1. 

In the limit of large J' the equivalence of the t-J ladder 
with a Luther-Emery liquid can clearly be seen. Going 
to isotropic coupling the spin gap remains finite and the 
only low-lying excitation is the collective charge mode, 
as we will show in the following sections. Thus also 
at isotropic coupling the t -J ladder is still a Luther- 
Emery liquid. In Sec. VI B, we will develop another ap- 



proach which relates the long-range correlations to ther- 
modynamic quantities for more general J's, based on a 
bosonization of density fluctuations. 



C. Phase separation 

Finally we study the occurrence of phase separation 
at isotropic coupling J — J' . We estimate the onset of 
phase separation by determining the coupling J at which 
the compressibility k diverges. The compressibility per 
site can be calculated as usual 



P 



2 ^Mp) 

dp^ 



(10) 



where e(p) is the energy density per site of the ladder 
with a particle density per site p = N/{2L). 
In a finite system usually the discrete version 



N 
2L 



2 r 



E{N + 2;L) + E{N -2;L)- 2E{N; L) 



ill) 



is used, where E[N; L) is the ground state energy of the 
finite system with N particles on the ladder with L rungs 
(Volume 2L). At small hole doping however this proce- 
dure may not be reliable due to finite size effects caused 
by frustration on small lattices. To see this let us con- 
sider the 1/ = 8 ladder doped with zero, two or four holes. 
In the undoped case there are 8 spins on each leg of the 



ladder. Two holes doped into the ladder will predomi- 
nantly go onto different legs and there will be seven spins 
per leg. Thus the antiferromagnetic configuration on the 
legs will be frustrated. For four holes there will be six 
holes on each leg and the system is again not frustrated. 
Conversely on an L = 9 ladder the undoped ladder and 
the ladder doped with four holes will be frustrated, while 
the ladder doped with two holes will be non-frustrated. 

We have thus used a different formula to calculate the 
compressibility at small hole doping. We calculate the 
ground state energies for an L = 8 ladder doped with 
A^h — and 4 holes and for an L = 9 ladder doped with 2 
and 6 holes. In all these cases the ladder is not frustrated. 
Then we estimate the compressibility from these energies 
using finite differences similar to the above Eq. (|ll|). 
In the thermodynamic limit L oo both formulas give 
the same result, as the frustration appears only on small 
lattices. 

While the finite size effects are quite small at low elec- 
tron densities they are much larger at small hole dopings 
due to frustration mentioned before. The estimated er- 
rors on the phase separation line may thus be much larger 
there, about ±0.2^. 

A comparison with the results obtained with open 
boundary conditions (OBC) confirms our results. Only 
at small doping the OBC results are not reliable since 
there the holes are trapped on the ends of the chain. 

Figure shows the phase separation line for the t-J 
ladder for J = J', in the. J-p plane. Note that, opposite 
to the single chain case,EZI the onset of phase separation at 
small hole doping is at lower values of J/ 1 than at small 
electron concenteations. This resembles the behavior in 
two dimensions ,B3 although the precise position of the 
phase separation line in two dimensions has not yet been 
established. 



IV. PROPERTIES OF A SINGLE HOLE 

In the previous section we have discussed the ground 
state of the ladder doped with two holes. In order to 
understand the low energy excitations of the ladder it is 
useful to study the one-hole problem first. 

As mentioned above the limit J' ^ J, t is a good start- 
ing point to explore the t-J ladder. There are nine dif- 
ferent states, depicted in Fig. ||. A single electron goes 
either into the bonding or antibonding orbital 
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V2 



1 / 



V2 



(c^ - ) 



(12) 



with energy =pi', respectively. Two electrons on the rung 
are either in the singlet state with energy — J' or in one of 
the three triplet states with energy 0. The singlet state 
expressed in bonding and antibonding orbitals is 
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— (c^ - - — (b^ - aM 



V2 



(13) 



Similarly the three triplets can be expressed as combina- 
tions of one bonding and one antibonding electron; 



atfet, — = ( alb 
T T' ^2 V T 



1 



-a\b] 



(14) 



Figures y(a) and (b) show the one-hole spectra for L = 
8 for large inter-chain coupling J' / J — 10, calculated 
by exact diagonalization for J/t = 0.3, J'/t — 3 and 
J/t = 0.5, J'/t = 5 respectively. 

A hole on a single rung can be either in the bonding 
or the antibonding orbital. One hole doped into the half 
filled ladder will thus be either in a bonding or antibond- 
ing state, depending on the parity symmetry of the total 
ladder [see Fig. ^(b)] . This hole can propagate along the 
ladder with a hopping matrix element t = +t/2 in first 
order perturbation theory. Thus the low energy states 
are two bands of holes in the bonding and anti-bonding 
orbitals. They are split by the energy difference 2t' of 
the bonding and antibonding states. These two bands 
can clearly be seen in the spectra [Figs, ^(a) and (b)]. 
The minimum of the bands is at fc^: = tt, since the hop- 
ping matrix element for holes i > 0. The bandwidth of 
both bands is 4t = 2t in the limit J' ^ J, t. At finite J 
the bandwidth is reduced due to hybridization with the 
higher excited states. 

Decreasing J' to the isotropic points J = J' changes 
the dispersion of these bands [see Figs.|(c) and (d)]. At 
low energies we can still see the bands of holes in the 
bonding and antibonding orbitals. These bands evolve 
continuously from the large J' limit. The minima of the 
energy bands are not at = or = but at a 
large momentum « (±|7r, O) for the bonding and 
k;^ ~ (^tfTT, tt) for the antibonding band. We can fit the 
low-lying hole bands to a dispersion of the form 



E{kx) = Eo + AE + ai cos k^ 

+a2 cos 2fcj. -I- aa cos Sk^, 



(15) 



corresponding to nearest neighbor (ai), next-nearest 
neighbor (q!2) and third-nearest neighbor (as) hopping. 
Eq is the ground state energy of the undoped ladder and 
AE the shift in energy of the center of the band upon 
doping. In Fig. |o| we show the bands and the excellent 
fit. The parameters are shown in Table |. 

The changes in the hole dispersion with decreasing J' 
are summarized as follows: 

(i) The center of the bands shifts downwards by AE < 
0, compared to the undoped ladder. The energy gain for 
one hole in the case of J = would be just the kinetic 
energy —t'. When J > we lose magnetic energy by 
introducing the hole. The energy gain is therefore smaller 
at larger J/t, as we can also see from the fit parameters. 

(ii) The hole bands are narrowed compared to the large 
J' limit. In that limit the bandwidth of the hole bands 



was 2t. This bandwidth is renormalized by the stronger 
polarization effects at isotropic coupling, and it is now of 
the same order as the magnetic energy J, instead of the 
kinetic energy 2t. 

(iii) The dispersion changes as longer range hopping 
processes (a2, as) are introduced with decreasing J', and 
the minima move away from fc^, = tt. The minima of 
both bands are very close in energy, again in contrast 
to the stro ng co upling region where they are split by 
2t' . In Sec. VII we will identify the minima kp and kp 
with the Fermi points of the bonding and antibonding 
quasiparticle bands. 

Another interesting question is the behavior of the free 
spin that is left over after one hole has been doped into 
the ladder. In the t-J chain the spin and charge exci- 
tations are carried by different soliton excitations which 
are far apart in space from each other. This is a typical 
feature of spin-charge separation and such a system is 
called a Luttinger liquid. In a Fermi liquid on the other 
hand they are bound and the excitations are described 
by quasiparticles carrying both charge and spin. 

We have calculated the hole-spin correlations to an- 
swer the question if spin-charge separation occurs in the 
ladder. The real space correlations 



(n„,a(j+r)5:,(j)). 



(16) 



are shown in Fig. This correlation function is nonzero 
for the ground state in the subspace of S'foj = ^ since 
there remains one spin unpaired. 

The result shows that the hole is tightly bound to 
the remaining free spin. At strong inter-chain coupling 
J' 3> J, i it is again predominantly on the same rung. 
At isotropic coupling the spin-hole bound state is more 
extended. These spin-hole bound states thus carry both 
charge and spin. In this sense they are similar to the 
quasiparticles in a Fermi liquid. This is in contrast to the 
spin-charge separation in the single chain. We will there- 
fore call the single holes bound to the free spin "quasi- 
particles" , although the system has a spin gap. 



V. EXCITATION SPECTRA OF THE LADDER 
WITH TWO HOLES 

A. Excitation Spectra 

The ground state of the ladder doped with two holes 
is, as discussed above, a bound state of the two holes. 
This bound pair coherently propagates along the ladder, 
giving rise to the lowest-lying band. When J' ^ J, this 
band, spin-singlet charge excitations, is clearly seen in 
the numerical results as shown in Fig. |l^. 

The higher energy excitations are again understood 
simply in the large J' limit. An essential difference from 
the lowest-lying singlet band is that two holes are now 
separate rather than forming a bound pair. Being sep- 
arate, they can gain a larger kinetic energy, but only in 
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return for a even larger cost of exchange energy J' as 
one more singlet rung is broken. Thus there are continua 
of scattering states of the two holes ( "quasiparticles" ) at 
higher energies. Since the residual interactions between 
the two quasiparticles are weak, the energy is almost de- 
generate between the S = and S = 1 spin subspaces. 
On the finite lattice we naturally do not see a contin- 
uum of scattering states, but only several discrete bands. 
These bands, and the fact that the energies of the triplet 
and singlet are nearly degenerate (up to boundary effects) 
can be seen in the spectra. 

There are various combinations of the two quasipar- 
ticle bands in the two quasiparticle continuum of states. 
The lowest are scattering states of two bonding quasipar- 
ticles, with ky — 0. Higher states are scattering states 
of one bonding and one anti bonding quasiparticle. Hav- 
ing the same transverse momentum, ky = t: these states 
hybridize with the one-magnon excitations in the spin 
background, resulting in a "bound state" below the two- 
particle continuum. This is nearly dispersionless and 
clearly seen in the spectrum. The continua of two an- 
tibonding quasiparticles is much higher in energy and 
not included in the figure. 

At higher energies there are spin excitations in the spin 
background. They are not described by the quasiparticles 
and will be discussed in more detail in the next section. 

The results in the isotropic case, J'/ J = 1, is shown in 
Fig. The energy spectrum is more complicated but 
the above description still holds qualitatively. 

The ground state is still the bound hole pair. It moves 
coherently along the ladder, yielding a gapless band of 
singlet charge excitations. The band has a linear disper- 
sion around k = (0, 0) compared with quadratic in the 
large J' case. An important point is that despite its com- 
plicated dispersion the low energy part is well separate 
from the other excitations. Therefore also at isotropic 
coupling the only low-energy excitation is the collective 
charge excitations, and we may identify the isotropic t- 
J ladder as Luther-Emery liquid. We will discuss the 
essential role of these c harge fluctuations concerning su- 
perconductivity in Sec. VI B| . 

In addition to the gapless band of charge fluctuations, 
there are various local minima at higher energies. How- 
ever, they can be explained by taking account of the non- 
monotonic dispersion of the one-hole spectra shown in 
Fig. ^ and our quasiparticle picture still holds. More 
specifically, there are four local minima in the single- 
hole spectra at k « (±|7r,0) and k « (±|7r,7r), which 
have nearly the same energy. Thus to construct low- 
energy two-quasiparticle excitations there are many pos- 
sible combinations of different minima of one-particle 
states as discussed below. This is the origin of many 
local minima in the two-particle spectra and this is con- 
firmed by the fact that the dependence on the boundary 
condition is consistent with this picture. 

The ground state of the two-hole spectrum can be con- 
structed from holes near zLkp in the single-hole spectra. 
The other minima can be explained similarly. The mini- 



mum in the PBC spectrum at kx = ^tt can be identified 
with the 2kp excitation, where a particle moves from one 
Fermi point to the opposite one. (2x^Tr= ^tt). When 
using APBC we do not have the fc-value of ±f tt in the 
single-hole spectrum on an L = 10 ladder. The closest 
/Caj-points are ±^7r and ±^7r, leading to minima at ±7r 
and ±|7r in the two- hole spectra. Another feature that 
can be explained from the single hole spectra is the min- 
imum at k = (tTjTt) (odd parity, k^ — tt). Using PBC 
this state is obtained with one hole with k = (|7r, 0) and 
one with k — (|7r,7r), leading to the minimum at (tTjTt). 

Using APBC's we can combine one hole at k = (j^t^, i") 
with one at either k = (^tt, tt) or k = (^7r,7r). As 
mentioned above these two states are higher in energy 
than the minimum at ifTi" and very similar in energy. 
Therefore we expect two states at k = (|7r,7r) and k — 
(tt, tt) which are similar in energy but at a higher energy 
than the corresponding states with PBC. This is exactly 
what we observe. The odd parity states near fc^; = can 
be constructed similarly. 

The qualitative features of the three-hole excitation 
spectrum can again be explained similarly. 



B. Spin Excitations 

One of the most interesting properties of the t-J lad- 
der is that there are two distinct types of spin excita- 
tion. Although it is most easily seen in the large J' limit, 
the qualitative distinction remains down to the isotropic 
point. 

The first type is the collective magnon excitations 
inherited from the undoped spin ladder. One of the 
electron- filled rungs is now excited to a spin triplet. This 
local excitation is what we call "magnon" and it propa- 
gates coherently along the ladder, leading to an energy 
dispersion with respect to kx- For a detailed investiga- 
tion, we have examined the two-hole spectra in more de- 
tail and have calculated the spin-spin and spin-hole cor- 
relations of the low- lying triplet states. We find that the 
magnon excitations of the Heisenberg ladder evolve con- 
tinuously upon doping. However there discontinuously 
appears a new kind of spin-triplet excitation at lower en- 
ergies, which is not present in the undoped ladder. 

The lowest excitation is a different type for which 
quasiparticles play an essential role. Therefore, the spin 
gap, defined as the excitation energy to the lowest triplet, 
is a discontinuous function of the hole doping at S — 0. 
This new type of excitations consists of breaking a pair 
of holes into two separate_quasiparticles, each carrying 
charge + \e\ and spin 1/2. When the two quasiparticles 
are both in the bonding orbital, their lowest energy, at 
k = (0,0), is lower than the lowest magnon excitation. 
The additional energy gain is easily understood, since the 
two separate holes have a larger kinetic energy of the or- 
der of t while the gain of the magnon kinetic energy is 
the order of t^/J'. As was shown in Fig. ||, down to the 
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isotropic point the lowest two-quasiparticle excitation is 
lower in energy than the lowest magnon excitation. 

The above picture is confirmed by comparing the cor- 
relation function between the two different states. The 
dynamical spin structure factor gives another confirma- 
tion. 

Figure |lj shows the equal-time correlations of the two 
holes, (nh(r)nh(O)), and of spin and hole, (5^(r)nh(0)), 
calculated for these two types of spin-triplet excitations. 
The latter quantity is nonzero since the states with 
= 1 are used in the calculation. When J' /J = 10, 
the two holes are separate in space in the lowest state, 
while tightly bound in the other state. The position of 
magnetic excitation is, on the other hand, close to the 
hole position in the lowest state, while they are far apart 
from each other in the other state. These two behaviors 
are what is predicted by our picture explained above, and 
despite modification in small detail they are qualitatively 
consistent even at the isotropic point. 

Our numerical results confirm that the lowest triplet 
state is the quasiparticle excitation where the bound hole 
pair breaks up. The two holes repel each other and the 
hole-hole correlations have the maximum at the largest 
distance L/2. Each of the holes is bound to a spin- 1/2, 
as can be seen from the hole-spin correlation function. 
This state is sketched in Fig. ^(d). 

As a typical magnon excitation, we show the correla- 
tions for the state at k = (vr, tt) that has the main spectral 
weight in the dynamical spin structure factor which will 
be discussed soon. The hole-hole correlations are simi- 
lar to the ground state. The hole-spin correlations show 
that the triplet carrying spin current is far away from 
the hole P^r. This state is shown in Fig. |^(e). Mean 
field theoryO predicts only this magnon excitation, which 
evolves continuously from the Heisenberg ladder. 

In neutron scattering experiments the relevant quan- 
tity is the dynamical spin structure factor 



>S(q,wo) 



S^\G.S.) diEr^-Ecs.-LOo), (17) 



where \n) is the complete set of eigenstates with energy 
En, |G.S.) is the ground state with energy Eq.s. and 



1 



r QZ 



(18) 



Figure shows 5(q, w) calculated for the Heisenberg 
ladder using the Lanczos diagonalization combined with 
the continued fraction method.t3 It can clearly be seen 
that the dominant contributions arise from the collective 
excitations near q = (tt, tt) and there is very little weight 
in the continuum of spin excitations at higher energy. 

In the doped case the two types of spin excitations have 
different contributions to the dynamical spin structure 
factor. It can be seen in Fig. |6| that the continuum of 
spin excitations move towards lower energies. Most of the 
weight is in the magnon excitations, consistent with the 



mean-field theoryjj and there is very little weight in the 
lowest triplet excitation consisting of the two separate 
quasiparticles. 

As in the large J' limit the lowest triplet excitation 
with Qy = t: is a bound state of a spin triplet and the 
hole pair. At q = (tt, tt) this state has no spectral weight, 
while most of the weight is in the second excited state, 
which has the triplet separated from the bound pair. This 
is a finite size effect of the two-hole system. The reason 
is that at = or TT we have an additional symme- 
try, reflection invariance in the ladder direction. The 
parity under these reflections is different for the ground 
state and the triplet-hole pair ground state, leading to 
the vanishing weight. 

A significant difference between the two types of spin 
excitations is that the largest number of "quasiparticle" 
excitations is limited by the number of holes. We need at 
least two holes to create such an excitation. The number 
of possible excitations is thus proportional to the hole 
doping S. On the other hand, the magnon excitations 
can be excited at any rung where there are no holes. 
The number of these excitations is thus proportional to 
1 — 6 instead, much larger for a small doping S. Therefore 
with decreasing temperature the susceptibility will show 
a large exponential drop at temperatures of the order of 
the gap of the undoped system (T ~ 0.5 J), followed by 
a small drop at temperatures around the spin gap of the 
doped system. 

To summarize we can describe the spin excitations of 
the t-J ladders by a simple picture: quasiparticles mov- 
ing in a spin liquid background. In the ground state 
the quasiparticles are paired. In the excitation spectrum 
two types of excitations can be distinguished. The first 
corresponds to the breaking of a pair of quasiparticles. 
This excitation has the lowest energy, but its number 
is limited by the number of holes. Of more importance 
for measurements of the susceptibility or inelastic neu- 
tron scattering experiments is the second type, which are 
magnon excitations in the spin liquid background. They 
evolve continuously from the undoped Heisenberg ladder. 
Although the gap for this type of excitation is larger it is 
more important since we can excite more of these exci- 
tations. Also the weight in the dynamical spin structure 
factor is larger. 



C. Charge Excitations 

Similarly to the calculation of the dynamical spin 
structure factor we calculate the dynamical charge struc- 
ture factor defined by 

AA(q,a;o) = ^| (n|pq|G.S.) |',5(^;„ - ^;g.s. - ^o), (19) 

n 

where pq is the Fourier transform of the local density 
fluctuation around the average density p. 
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The result of our calculations is shown in Fig. |l^. The 
main contribution arises from the coherent motion of the 
hole pairs. This leads to large peaks at low energies near 
q = (0, 0). The rest of the weight is distributed incoher- 
ently over a wide region at rather high energies of order 
At. It arises from interactions of a single hole with the 
surrounding spin background-. This is similar to results 
obtained for two dimensions. lj 

Recently Tohyama et al. have calculatpi 5(q, uj) and 
A^(q, uj) for single chains and 2D clusters.LJ In the single 
chain they find that, as expected from Luttingcr liquid 
theory, the charge and spin excitations are decoupled in 
the low energy region. The dynamical charge structure 
factor is very similar to that of spinless fermions, con- 
sisting of large peaks at the energies expected from the 
cosine band of the spinless fermions. 

In two dimensions they find different behavior. In the 
spin structure factor nearly all of the weight is in a few 
sharp peaks at low energies. In the charge structure on 
the other hand the main weight is at relatively large ener- 
gies, of the order of several t, and it strongly broadened. 
This is indicative of strong spin-charge interactions. Only 
at some special q-points are there peaks at low energies. 
The dynamical charge structure factor, A/'(q, w) for the 
ladder shown in Fig. |l^ differs a lot from the single chain 
but resembles the 2D cluster. At large Qx, we see the 
peaks at large energies ~ 2t-At. At small Qx in the Qy = 
sector, A/'(q, uj) for the ladder is dominated by the col- 
lective mode of the hole pairs, and its behavior differs 
substantially from spinless fermions. Another similarity 
to the 2D system is found in the dynamical spin structure 
5(q, Lu) where the main weight is in peaks with energies 
~ J (see Fig. |l6|). We conclude from this comparison 
that the 2D clusters and the ladders are closely related - 
a fact which points towards a "d-wave" paired state for 
the 2D clusters. A more difficult question is to what ex- 
tent this behavior of the finite 2D cluster is a consequence 
of the strong tendency of finite clusters to favor singlet 
ground states and to what extent it is representative of 
an infinite plane. 



duce the operator creating a singlet pair of electrons on 
sites r and r -I- d. 



1 / 



Pi. = ^ c' c. 



,t J 



-old 



(20) 



Using this definition of the pair operator we can calculate 
the superconducting order parameter 



= (G.S.,7Vh-2 holes] ^^Pjd 



Xd 



and its Fourier transform 



G.S., Nh holesj 



(21) 



Xk= 



^Xd' 



,jkd 



= (G.S., 7Vh-2 holes|P^__i^|G.S., A^h holes) (22) 

where Pk,-k is the Fourier transform of the real space 
pair operator with zero total momentum. 

Figure |l^ shows the superconducting order parameter 
for the L — 8 ladder with PBC's for several values of 
J/t. The most obvious properties are that the sign is 
opposite for pairs with ky = and ky = n and that 
the absolute values are very small near (0,0) and (7r,7r). 
This is similar to rf-wave pairing in a fully 2D system. 
The absolute value is largest near wave vectors which we 
have identified with the Fermi points of the quasiparticles 
k = (37r/4,0) « kf and k = (7r/2,7r) w k^. At large 
J/t the order parameter is very similar to the simple 
cos kx -— cos ky structure of nearest neighbor c?-wave pairs. 

Consider now pairs with a form factor F: 



Pi 



d 



F*P. 



r,d' 



where F is normalized 



(23) 



(24) 



VI. LONG RANGE CORRELATIONS 



A simple calculation shows that the superconducting "or- 
der parameter" 



A. Superconducting Correlations 



Xf^^Fa Xd 



(25) 



A highly intriguing and much debated subject of the 
high-Tc superconductors is the internal symmetry of the 
order parameter. In this section we study the internal 
structure of the pairs in the doped ladders. For the t-J 
model usually only nearest-neighbor pairs have been con- 
sidered except for a few cases. It is reasonable to assume 
that they are the dominant pair correlations. However 
more quantitatively we have chosen the optimal form for 
the pairs. 

Using the Lanczos algorithm we calculate the pairing 
correlation functions for different pairings. Let us intro- 



is maximized by the pair with the form factor 

Xd 



-Fd 



Xr 



1/2 ' 



(26) 



which is proportional to the order parameter xa- This 
optimal order parameter then is 



[El 



x= I2^lxdi 

d 



1/2 



(27) 
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In Fig. ^ we show this maximal order parameter for 
a L = 8 ladder as a function of J/t. A clear increase 
as a function of J/ 1 can be seen, similar to single chains 
and 2D planes where the superconducting correlations 
are also enhanced at larger J values. We note that the 
expectation value of the optimal order parameter is larger 
for four doped holes than for two holes. This is a con- 
sequence of the strong correlation effect, as any pairing 
order parameter or fluctuations should vanish as (5 — > 0. 
The mean-field theory including the strong correlation 
effects predicts a (5- linear dependence of the order param- 
eter as i5 ^ OH i.e., the superfluid density is proportional 
to the hole doping. Our results show that the value for 
four doped holes is less than twice the value for two holes. 
This may be because that in the four hole case the dop- 
ing, 5 = 0.25, is already large and out of the region of 
(5-linear dependence. We need further study using larger 
lattices to clarify this point more quantitatively. 



B. Mapping to a Luther-Emery Liquid 



are determined in this way, it is straightforward to cal- 
culate the power-law exponents of correlation functions 
as follows. 

The Efetov-Larkin Hamiltonian is actually identical to 
the single-component (i.e., spinless) Luttinger model di- 
agonalized b¥-,bosonization, which was studied by Mat- 
tis and Liebj23 Luther and Peschel,a and Haldanenj in 
details. The Hamiltonian (^) is immediately solved by 
rewriting each Fourier component with a boson operator. 



ns.k 



(29) 



then the result is 

Hel -5I^-I^I^I^fc + ?7K(^B"^B)'+«./^il- (30) 



In Sees. 1110, our numerical results have shown that 
all the spin excitations cost a finite energy and the only 
gapless excitations are charge fluctuations, i.e., coherent 
propagation of bound hole pairs along the ladder direc- 
tion. We mau-therefore say that the system is a Luther- 
Emery liquidcj from large J' limit down to the isotropic 
coupling. This observation indicates that the low-energy 
and long-wavelength properties of lightly doped t-J lad- 
ders would be entirely described by an effective contin- 
uum Hamiltonian in terms of charge degrees of freedom 
and this model correctly predicts long-range asymptotic 
behavior of correlation functions. The effective model is 
actually tha-hosonic Gaussian model proposed by Efetov 
and Larkin,L3 



Here Nb is the total number of bound hole pairs and 
vn — {t^K-q)^"^ is the charge velocity and vj = ttRbv^ is 
the current velocity associated with the number of 2fc|,^ 
excitations, Jb, obeying the universal relation, vnVj = 
v^. Here = nS = 7r(l— p) is the Fermi wave number of 
spinless fermions transmuted from the bosonic hole pair 
operators by the Jordan- Wigner transformation. The 
last two terms in Eq. ( ^0| ) describe non-bosonic excita- 
tions accompanied with the change in the quantum num- 
bers, A'b and Jb. Their importance was firai-pointed out 
by Haldane in his Luttinger liquid concept Hj 

The propagator of bound hole pair, i.e. pairing correla- 
tion, is then obtained as in the calculation of the Debye- 
Waller factor, 



G(r)= ns (e 



cx r 



'Hel= 



dx 



[nB{x) 



where the exponent is 



(28) 



Kp = 2t:KbVs = TT/d Vs k, 



(31) 



(32) 



where nB{x) represents the density of bound hole pairs 
at the rung x with the average value ub = Nb/L — 6, 
and Ob{x) is its conjugate phase obeying the canonical 
commutation relation, [nB{x),9B{x')] = i6{x — x'). As 
shown by the fact that the q = mode of the first term 
in Eq. (pa) denotes the change in the ground state en- 
ergy associated with the number of bound hole pairs, 
AE (X {1/2Kb){A.NbY , the parameter Kb is given by 
the compressibility of the hole pairs, which will be shown 
afterwards. The second term, on the other hand, de- 
scribes the energy change associated with current, related 
with the sound velocity Vs- Instead of a direct calcula- 
tion of the dispersion relation, the sound velocity can 
also be obtained numerically by applying an Aharanov- 
Bohm flux penetrating in the center of the ladder with 
PBC's, since the flux induces a finite current along the 
chain direction. Once the two parameters, Vs and Kb, 



where p is the electron density per site and n is the 
electron compressibility. The relation for the ladder, 
Kb — K/9^/2, is used to obtain the second equality. This 
result differs from the one for the chain, which has the 
numerical factor 7r/2 instead. It is important to empha- 
size that this result is the Luttinger liquid parameter, 
Kp = 2{vj /vj^y^^ , for the single-component boson sys- 
tems. 

Density-density correlations are characterized by a 
power-law decay at the wave number 2k^ . Since is 
twice the average of the Fermi wave number of the orig- 
inal electron bonding and antibonding bands, kp = ^p, 
the wave number 2k^ is actually Akp in the original pic- 
ture. The exponent of the 2kfF=Akp CDW fiuctuations 
is calculated using a_density-phase duality as was done in 
the original paper.t3 The density operator has a short- 
range 2kf CDW order, but since the value of 2k% is 
locally determined by the density at each position, it is 
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fluctuating around the average value 2tt5 associated with 
the density fluctuations. We therefore may write as 



nsix)— ub cos[2fc|?"x + 26j{x)] , 



d9j _ 



(33a) 
(33b) 



Rewriting the Hamiltonian (28) with Oj and its conjugate 
operator nj = ^V^s, we again obtain a Gaussian model, 



dx I vj[nj{x)]'^ + VN 



VOjix) 



(34) 



where [nj(a;), 0,/(a;')] — iS{x — x'). This duality of the 
charge and current operators was also emphasized in the 
Haldane's paper.Cj Since the density-density correlation 
is written as an exponential of 9j in this representation, 
its calculation can be done similarly as before, 

NB{r)= (ns(r)nB(O)) «nB2cos(2fcf r)(e2»[«-'W-«^(°)l) 
cx cos(2fc|M^~'^''- (35) 

Instead of calculating the sound velocity, the value of 
Kp is more accurately obtained through the Drude weight 
defined by 



D 



2 5(^2 



(36) 



= 



where EQ s {(j)) is the total energy of the ground state 
when the flux (jj is penetrating. Assuming the holes prop- 
agate in pair, the relation between D and for the lad- 
der is-Cjbtained and the result agrees with the one for the 
chainLJ aside from its numerical factor. 



D 



■V.J 



KpVs ■ 



(37) 



Combined this with Eq. (p2[), we finally obtain the for- 
mula of the correlation exponent in terms of the com- 
pressibility and the Drude weight: 



(38) 



It is noted that this expression is identical to the one for 
the chain in terms of k and D, while they differ by factor 
2 in terms of n and Vg. 

In this way, when the energy scale concerned is smaller 
than the spin gap, we can predict long range asymp- 
totic behavior of the correlation functions based on the 
Efetov-Larkin effective model, and calculate correlation 
exponents once the two parameters, Vs and (or equiv- 
alently k and D), are numerically determined. 

Figure ^ shows the correlation exponent calculated in 
this way for the isotropic ladder {J' / J = 1) with L = 7 
and two holes as a function of J. This corresponds to the 
electron density, p = | = 0.857. The necessary quantities 



are carefully calculated by using the Lanczos diagonal- 
ization. The compressibility, k, is determined from the 
ground state energy of unfrustrated systems in the sense 
explained before. The Drude weight, D, is then calcu- 
lated by Eq. (|36|). Combining these two, Kp is obtained 
via Eq. dH). 

Near the phase separation boundary, Kp grows rapidly. 
This is owing to the divergence of the compressibility at 
this boundary. In other words, the collective charge ex- 
citations become softening and the superconducting fluc- 
tuations are enhanced correspondingly. 

Recently Hayward et al.Ea directly calculated vari- 
ous correlation functions for the t-J ladders by using 
the density-matrix renormalization-group method. At 
J/t — J' /t — 1 and p = 0.8, they found a power-law 
decay of the pairing correlations with the exponent close 
to unity. This exponent corresponds to Kp ^ 1, which is. 
larger than our estimate at the same J and J', Kp « 0.7.0 
The discrepancy may be owing to the larger electron den- 
sity in our calculation, p = | = 0.857, but it is not clear 
if this suffices to account for the difference until new cal- 
culation is carried out at the same density. Hayward et 
al. alsOpCeport data on the density-density correlation 
functionE^ but the exponent of the expected power-law 
decay of the oscillatory term cannot be easily extracted. 



VII. ONE-PARTICLE EXCITATIONS 

Finally we discuss the one-particle Green's function 
where we can see the quasiparticle excitations directly. 
The electron and hole parts of its spectral function are 
deflned as 

A,,AK ^)= Y.\^n, 2L-l|c[jG.S., 2L-2) |' 

n 

xd{Lj - S„(2L-1) + ^G.s.(2L-2) + /^), 
Ah,,(k, ij)= ^1 (n, 2L-3|ck.|G.S., 2i-2) |' 

n 

x5{uj + En{2L-3) - Eg.s.{2L-2) + p), (39) 

where \n,N) is an eigenstate for N electrons with the 
energy En[N) and |G.S.,7V) denotes the ground state 
with N electrons. Positive (negative) energies correspond 
to the electron (hole) part. The chemical potential is 
defined by 



i[i;G.s.(2i-l)-£;G.s.(2L-3)]. 



(40) 



The results for L = 10 and the isotropic couplings, 
J'/J = 1, are shown in Figs. |l|for J/t=0.3 and 0.5. The 
wave vectors along the ladder kx = {n: integer) are 
for PBC's and kx = ^{n+\) for APBC's. The ground 
state energy Eq^{2L — 2) and the chemical potential p. in 
Eq. ( p9| ) are the average over both boundary conditions. 
The results for L — lOr-are very similar to our L = % 
results published before. B 
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The spectral function has large weights for the bonding 
(B) {ky = 0) and antibonding (A) {ky = n) orbitals only 
near the Fermi energy w = 0, and they seem to constitute 
quasiparticle bands. Away from the Fermi energy, the 
individual quasiparticle peaks arc much less prominent 
and there is an incoherent part with an energy of the 
order of t. 

The quasiparticle part of the spectrum is consistent 
with ±he mean-field theory based on the d-wave RVB 
state.S The undoped ladder consists of local singlets on 
the rungs. Such a singlet is the superposition of two 
electrons in the bonding orbital and two electrons in the 
antibonding orbital, — a|aj |0). Holes doped into 
the half-filled ladder will go predominantly into the anti- 
bonding orbitals to gain a larger kinetic energy along 
the rung direction. The bonding band is occupied by 
more electrons, while the antibonding band is occupied 
by less electrons. The quasiparticle with energy clos- 
est to w = has a wave vector nearest to the original 
Fermi wave number, kp- {kx — ^ for bonding and 
kx ~ ^ for antibonding). Because of the band split- 
ting, fc^ > kp, but the Luttinger sum rule is satisfied, 
kf+k^ = {1-S)tt. This means the Fermi volume is large, 
proportional to the electron number rather than the hole 
number, and this is consistent with pJiipto-emission exper- 
iments on cuprate superconductors.EZI It is important to 
notice that the quasiparticle peaks near the Fermi energy 
have their counterparts on the opposite side of the Fermi 
energy. An electronic quasiparticle peak at energy uj > 
has a shadow hole peak at energy around —lo < 0, and 
vice versa. These peaks indicate that the quasiparticle 
excitations are those of the Bogoliubov quasiparticles as 
in BCS theory, i.e., mixture of an electron and a hole 
(aj^ — ""kCkj + ■f^kC-kx)- The weights in the electron and 
hole parts are proportional to l^kP and |wkP- They are 
hole-like around fc^; = and electron-like around k^ = tt 
for both the bonding and antibonding bands. 

There exists a finite energy gap in the quasiparticle 
spectra. The electron and hole branches both come close 
to the Fermi energy ai k^ ~ -f, but instead of pass- 
ing through they move away from it. The energy gap 
for J/t = 0.3 (0.5) is O.lSt (0.29t) at k = (^,0) and 
0.22^ (0.39i) at k = (^,7r). This corresponds to a 
quasiparticle gap 2Aqp ~ 0.13t ~ ^ for J/t — 0.3 and 
2Aqp ~ 0.29i ~ ^ for J/t = 0.5. 

It is interestiiW-|to note that the calculations of A(k, uj) 
in 2D clustergL3^t2l show similar behavior for k-points not 
along (1,1) but no shadow peaks for k || (1,1), indicating 
d2,2_j^2 -pairing also. 

Figure |2^ shows the spectral function of the one- 
particle Green's function of two holes in an L = 10 ladder 
summed over all wave vectors: Ao,ct(w) — J2k ^c,o-(k, uj), 
and Ah,cr(w) — Ah,cr(k, w). This quantity is the lo- 
cal density of states to add and remove electrons. In a 
strongly correlated system, the sum rules on the weight 
(i.e., the integrated values) of Ac{uj) and Ah(w) are very 



different since the former is given by the number of empty 
sites and the latter by the number of filled sites (or equiv- 
alently the number of holes and electrons, respectively). 
In a Fermi liquid the values of Ac{uj) and Ah(w) for small 
values of jcij — /i| are determined by quasiparticle weight 
at the Fermi energy and are continuous. It is interesting 
therefore to note that Fig. |2^ shows approximately simi- 
lar values for Ae{uj — n) and Ah(/i — w) around lu ^ fi, but 
the sum rule on the total weight is satisfied through the 
large weight in incoherent excitations in Ah(a') at higher 
energies, |a; — /i| > J. The strong correlation condition 
is reflected in the much smaller total weight in Ac{uj) 
which comes about through an effective cut-off in energy 
on Ac{uj). In this respect the system in energy space is 
similar to a lightly hole doped band insulator although 
as we discussed earlier the location in k-space of the co- 
herent quasiparticle peaks corresponds to a large Fermi 
surface to add and remove electrons. 

The momentum distribution for electrons, n'^{^) = 
('^k.ffCk.A and for holes, n^{k) = (ck,crCk_^), is shown 
in Fig. Note because of the strong correlation con- 
dition, these do not add to one but instead their sum is 
given by n° (k) -I- n^(k) — ^{1 + S). The strong correla- 
tion condition is also evident in the reduced magnitude 
of the variation of n°(k) as a function of k. Nonethe- 
less the presence of an apparent "Fermi surface" in the 
center of the Brillouin zone is clear, consistent with the 
dispersion relations of the coherent quasiparticles. The 
difference between the bonding (ky = 0) and antibond- 
ing (ky — tt) bands and the reduced occupation of the 
antibonding band arise from the energy gain in placing 
the doped holes preferentially in the antibonding band. 

VIII. CONCLUSIONS 

The results of our Lanczos diagonalizations confirm 
earlier studies which concluded that lightly doped two- 
leg ladders belong to a different universality class from 
single chains. The latter are Tomonaga-Luttinger liquids 
with gapless and separated spin and charge excitations. 
The ladder in contrast has a finite gap in the spin excita- 
tion spectrum and gapless excitations only in the charge 
sector. The low energy excitations evolve continuously 
from the limit of strong inter-chain exchange coupling 
(J' ^ J, t) and the simplicity of that limit allows a clear 
interpretation of our results. 

At large J' the dispersion relation of a single doped 
hole consists of two cosine bands corresponding to bond- 
ing and antibonding states on a rung. Lowering J' to 
the isotropic limit {J' / J = 1) and setting both J,J'<t 
changes the dispersion relation substantially. The coher- 
ent parts of both bands are centered at energies ^ — 1.5i 
but the width is ~ J only. The spin and charge com- 
ponents are still bound but more loosely and the large 
magnetic polarizability of the spin background intro- 
duces longer-range hoppings. Remarkably the form of 
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the bands resembles the noninteracting band structure 
so that a photoemission experiment which removes elec- 
trons would measure in effect a large "Fermi surface" 
with bonding and antibonding pieces. In this regime the 
quasiparticle propagation is strongly influenced by the 
coupling to magnetic excitations. 

When two holes are added they bind together on a 
single rung at large J' and remain bound although the 
size of the bound hole pair increases as J' approaches 
J t/2. Moreover the qualitative features of the density- 
density A/'(q, w) and spin-spin structure factor iS(q, a;), 
which are easy to understand at large J', remain simi- 
lar as J' approaches J ^ t/2. A/'(q, w) near q = (0,0) 
is dominated by the low energy mode associated with 
the motion of hole pairs. At large qx, J^{q,uj) has a 
broad peak at high energies {'^ At) similar tOj-that found 
by Ohta, Eder and Maekawa for 2D clusters.Lj They in- 
terpreted this as local excitations of single holes in the 
magnetic cloud or spin bag.Lj The dynamical spin struc- 
ture factor, 5(q, w), also resembles 2D clusters and not 
ID chains when we compare to the results of Tohyama, 
Horsch and Maekawa.Ell The major weight is at energies 
~ J. The spin gap evolves discontinuously upon dop- 
ing through the new quasiparticle excitations that can 
be made by breaking a hole pair into two separate single 
holes. However the major weight of the spin excitations 
remains in the collective magnon mode whose dispersion 
evolves continuously from the 6 = limit, although it is 
influenced by the continuum of quasiparticle excitations. 

We have also investigated the one-particle spectral 
functions to add and remove electrons from the two-hole 
ground state. These show clearly the unusual nature 
of this "Fermi liquid". When electrons are removed (or 
holes added), the spectral weight is spread over a large 
energy region (^ 6t) , but the coherent part is limited only 
to energies ~ J below the Fermi energy fi. The energy 
dispersion relations show a large apparent Fermi surface 
for the coherent quasiparticles and which matches onto a 
similar one for adding electrons at energies greater than 
fi. These k-space features resemble a metal with a large 
Fermi surface. The property that resembles a lightly hole 
doped insulator is the energy dependence of the spectral 
weight to add an electron. This shows a low energy cut- 
off 6-6t) similar to a lightly hole doped band insulator. 
The result is an intriguing duality between metallic-like 
features in k-space and lightly hole doped insulating fea- 
tures in energy space. 

The overall properties of the lightly doped ladder place 
it in the Luther-Emery class rather than the Tomonaga- 
Luttinger class of ID systems. The low energy proper- 
ties of Luther-Emery liquids are described by interact- 
ing hard-core bosons as shown by Efetov and Larkin. 
In the present case the Efetov-Larkin bosons are bound 
hole pairs. Two features distinguish the t-J ladder from 
the usual Luther-Emery liquids arising from attractive 
interactions. One is the d-wave character of the pair- 
ing and the second is the presence of magnon excitations 
and limited quasiparticle excitations. Note the magnon 



excitations cannot be viewed as the collective mode of 
quasiparticles since the latter vanish as ^ ^ 0. The sys- 
tem is not a standard Fermi liquid, but rather is a new 
and interesting mixture of a dilute attractive Fermi gas 
in which the hole binding energy remains finite as (5 — > 0, 
and a dense Fermi liquid with an apparent large Fermi 
surface in k-space. 

Comparing the ladder with the results by Tohyama et 
al., we see that the ladder is very different from the single 
chain but similar to 2D clusters in many respects. Both 
in ladders and in 2D clusters d-wave pairs are found down 
to small J/t. The dynamical charge and spin structure 
factors look remarkably similar and the single-particle 
spectral functions indicate the existence of Bogoliubov 
quasiparticles with a finite superconducting gap. Thus 
we are lead to the conjecture that the t-J model on 2D 
clusters is a doped RVB spin liquid showing d-wave pair- 
ing, similar to the ladder. 
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TABLE I. Parameters for the fit of the lowest lying bands 
of the one-hole spectra to a dispersion of the form of Eq. (|l5[). 
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FIG. 1. The t-J ladder with two legs and L rungs. The 
couplings along the legs are t, J and those along the rungs t' , 
J'. 

FIG. 2. Ground state energies for the t ladder {t-J ladder 
with J = J' = 0) with two holes. Results are shown for sys- 
tems with L = 2, 3, . . . , 10 rungs and periodic (PBC) as well 
as antiperiodic (APBC) boundary conditions. The ferromag- 
netic state always has the lowest energy for L > 4 rungs. 

FIG. 3. Real-space spin correlations for the singlet ground 
state of the t ladder with two holes. L = 10 and PBC's are 
used. The two ferromagnetic domains can clearly be seen. 

FIG. 4. Graphical representation of the low-lying states of 
the t-J ladder in the strong coupling limit J' S> J, t. (a) 
The undoped case, (b) One hole goes into either the bonding 
orbital or the antibonding orbital on one rung, (c) In the 
ground state for two holes both holes are on the same rung, 
(d) Scattering states of two holes, (e) At higher energies there 
is the triplet excitation similar to the undoped ladder. 



FIG. 5. Binding energy of two holes, spin gap and en- 
ergy of the triplet excitation away from the bound hole pair. 
J/t = 0.3 and 0.3 < J' /t < 3.0. The size of the ladder is 
L = 8 rungs. 

FIG. 6. Hole-hole correlation functions for the ground state 
of the t-J ladder with two holes. J/t — 0.3, and J'/J=l and 
10. The size of the ladder is L = 10 rungs and APBC's are 
used, which have a lower ground state energy than PBC's. 
The inset shows the size of the bound hole pair, ^, in the 
two-hole ground state as a function of J' for different ladder 
sizes. 

FIG. 7. The line of phase separation in the t-J ladder de- 
termined from the coupling at which the compressibility di- 
verges. 



FIG. 8. The nine different states for a single rung. 



FIG. 11. The hole-spin correlations for the ground state 
of the ladder doped with one hole at J' /J = land 10. The 
ladder has L = 10 rungs and J/t = 0.3. The ground state has 
= 1/2. 



FIG. 12. Energy spectra for the t-J ladder doped with two 
holes for a large J' (J'/J = 10): J/t = 0.3 with (a) PBC's and 
(b) APBC's, and J/t = 0.5 with (c) PBC's and (d) APBC's. 
The size of the ladder is L = 8 rungs. The states are classified 
according to total spin S and parity. 



FIG. 13. Energy spectra for the t-J ladder doped with two 
holes at the isotropic point {J'/J = 1): J/t = 0.3 with (a) 
PBC's and (b) APBC's, and J/t = 0.5 with (c) PBC's and (d) 
APBC's. The size of the ladder is L = 10 rungs. The lines are 
only guides for the eye and do not necessarily connect related 
states. 

FIG. 14. Hole-hole [{a) and (c)] and spin-hole [(6) and 
(d)] correlation functions for the two triplet excitations. 
L = 8. Dashed lines are for the lowest triplet state and 
dashed-dotted lines for the lowest state with non-vanishing 
weight in 5(q= {■k,tt),lo). 

FIG. 15. Dynamical spin structure factor, 5(q, tj), for the 
undoped ladder with L = 10 rungs. Note that the scale is 
different for qy = n. 

FIG. 16. Dynamical spin structure factor, 5(q, oj), for the 
L = 10 ladder with two holes. J/t = J' /t = 0.3 (upper 
panels) and 0.5 (lower panels). PBC's are used. Note that 
the scale is different for around q = (tt, tt). 



FIG. 17. Dynamical charge structure factor, A/'(q, i^), for 
the L = ladder with two holes. J/t — J' /t — 0.3 (upper 
panels) and 0.5 (lower panels). PBC's are used. 



FIG. 18. Superconducting order parameter, Xk for Ni^ = 2 
calculated on the L = 8 ladder with PBC's for several values 
of J/t = J'/t. 



FIG. 9. Energy spectra for the t-J ladder doped with one 
hole. The case of large J' (J'/J = 10): (a) J/t = 0.3 and 

(b) J/t = 0.5 with L = 8. The isotropic case {J'/J = 1): 

(c) J/t = 0.3 and (d) J/t = 0.5 with L = 10. The results 
for kx = j^n with even n are for PBC's and with odd n for 
APBC's. 



FIG. 10. Fit of the lowest lying bands of the one-hole spec- 
tra to the form of Eq. (^. (a) J/t = J' /t = 0.3 and (b) 
J/t = J'/t = 0.5. The size is L = 10. 



FIG. 19. The order parameter for the optimal pair Eq. ( [2 7] ) 
as a function of J/t on an L = 8 ladder with both PBC and 
APBC. Shown are results for 2 and 4 holes. 

FIG. 20. Correlation exponent, Kp, for the L — 7 ladder 
with two holes. 
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FIG. 21. Spectral function of the one-particle Green's 
function, yl(k, w), for the L = 10 ladder with two holes, (a) 

J It = J' It = 0.3 and (b) J/i = J' It = 0.5. The width of each 
line represents the strength of the excitation. For w > we 
show the spectral function for adding one electron yle,cr(k, w), 
and for ci; < the spectral function for removing an electron 
Ah,CT(k,(i;). 

FIG. 22. Spectral function of the one-particle Green's func- 
tion for two holes summed over all wave vectors; Ah,cr(w) 
=Yl,]i^h,<y{^,i^) and Ac,cr{io)= A,<T(k,ci;). L = 10 and 
J It = J' It = 0.3. The oscillations at large \uj\ arc caused by 
nonconvergent Lanczos iterations at these energies. 

FIG. 23. The momentum distribution function for the 
ground state for 5 = 0.1 and L = 10. J/t = J' It = 0.3. The 
electron part, n%{\i)= {(^„c\s.^a), and the hole part, njj(k)= 
(ck.o-cj^ o-). The momenta, kx = 27rn/L, with integer n's are 
for PBC's, while those with half-integer n's are for APBC's. 
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